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1 Introduction 



Many of the Liouville integrable lattices, for instance the Toda lattice and the Heisenberg 
chain, remain integrable when imposing some boundary conditions other than the open or 
the periodic ones. Usually it corresponds to switching from the An-i classical root system 
to the other ones (-B„, C„, Dn, BCn, etc.) associated to an integrable lattice ^. It 
was an idea of Sklyanin ^ (see also |^]) to describe all such possible boundary conditions 
in terms of the representations of a new algebra 

R{u - w)r(^) {u)R{u + t;)r(2) (v) = r(2) {v)R{u + t;)r(^) {u)R{u - v) , (1.1) 

which was called afterwards the reflection equation algebra (see for instance ||5|). We 
prefer to give it another name the QISM II algebra (see §]), where QISM stands 

for the Quantum Inverse Scattering Method and II symbolizes the difference of it from the 
QISM I algebra, which is nothing but the familiar algebra given by the quadratic relation 

R{u - v)T^^^ {u)T^^'> (v) = r(2) (i;)t(^) {u)R{u - v) . (1.2) 

There are actually too many names for these algebras, and no one of them is becoming 
standard; that is why we are insisting on our own choices - the QISM I and II algebras. 
These names are short and quite characteristic too. We should also mention the work ||l2[ 
where a different version of the reflection equation (|1.1|) was considered. 



In terms of the integrable models the representations of the QISM I algebra (1.2) pro- 
vide us with integrable lattices (through the co-multiplication) while the representations 
of the QISM II algebra ( |1.1D describe the possible boundary conditions for such lattices. 
In this paper we use the new representation of the QISM II algebra found recently in |^ 
to obtain some new boundary terms for the known integrable lattices. These new models 
are generically formulated in terms of the XXX Heisenberg magnet interacting with two 
additional spins on each end. We consider as well the degenerate cases given by the cor- 
responding contraction procedure. We would like to stress that our results give for the 
moment the most general boundary terms for the Heisenberg magnet (expressed in terms 
of two additional spins on each end) and specialization will give already known boundary 
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conditions |T], ^ coming from the scalar or rank solutions of the QISM II algebra. De- 
scribing the situation in a bit more technical terms we can say that we are classifying all 
boundary conditions coming from the most general rank 1 representations of the QISM II 
algebra in the XXX case. All results in this paper are given for the Poisson algebra. The 
quantization is straightforward and will be published elsewhere. 



2 Rank 1 representations of the QISM II algebra 

In the classical limit the QISM II algebra (|1.1| ) becomes the following Poisson algebra 

{rW(u),r(2)(^;)} = [r{u-v),T^^\u)T^^\v)] (2.1) 

where T^^\u) = T{u) ® I, T^'^\v) = I T{v). Here denotes the usual tensor product 
and / is the 2x2 identity matrix. The T(u) is the 2x2 monodromy matrix depending on 
the complex spectral parameter u. The 4x4 matrix r{u) will in this paper be given by 

1 \ 
10 
10 

y 1 y 

This is the rational case, which corresponds to the XXX model |9|. 

Let T{u) be of the following form |^ 

A{u) B{u) 
C{u) D{u) 

with 



r{u) 



u 



T{u) 



(2.2) 



(2.3) 



Aiu) 
B{u) 



au + Aiu + ^0 H — ■ 



u 



D{u) = -A{-u) , 



pu^ + Bo , C(n) = + Co . 



(2.4) 
(2.5) 



Here a, /3, 7, and S are scalars, while Ai, Ao, Bo, and Co are generators of some algebra. 
The T{u) (p.3|)-( |2!5|) satisfies the following symmetry property 

T(-n) ~r-i(n). (2.6) 



Inserting this Ansatz for T{u) into the Poisson algebra (|2.l| ) leads to the following quadratic 
Poisson algebra A for the generators Ai, Aq, Bq, and Co 

{Ao,A^} = PC0-7B0, {Bo,Ao} = 2A^Bo-2p6, (2.7) 

{Bo,Ai} = 2aBo-2pAo, {Cq, Aq} = -2AiCo + 2j6 , (2.8) 
{Co,Ai} = 2jAo-2aCo, {Co, Bq} = AA^Aq - Aa6 . (2.9) 

The determinant of the monodromy matrix T{u) is the generating function for the center 
of the QISM II algebra (U). For the Ansatz (|1|)-(|2]|) 



det T{u) = -{a^ + p-f)u'^ + Qs^^ + Qo + ^ , (2.10) 



with 



Q2 = ^1 - 2aAo - /3Co - 7^0 , Qo = 2Mi - vl^ - BqCo . (2.11) 
Hence (2.11) gives two Casimir elements for the algebra (|2.7| )--( ^^ . 



The QISM II algebra ( p.lj) with the r-matrix ( |2.2[ ) admits the following scalar solution 

. a + d/u b 
K{u)=\ I, (2.12) 

c —a + d/u 

where a, h, c, and d are complex constants. This is just a special case of (|2.3D-( p^ ) with 
a = /3 = 7 = ^i=0. Then Aq, Bq, and Cq all commute. We may now combine the two 
solutions K{u) and T{u) of the QISM II algebra and define 

t{u) = tr K\-u)T{u) . (2.13) 



Note that Z : T{u) ^ T^{—u) is an automorphism of the algebra ( |2.1| ) with the r-matrix 
(^ . It is a property of the QISM II algebra [|] that 

{t{u),t{v)} = ^. (2.14) 



Hence the t{u) ( |2.13 ) is the generating function for the integrals of motion of an associated 



integrable system. For our special choices of K{u) and T{u) we get 



t{u) = {2aa + (3h + -ic)u^ + H - 2d^ , (2.15) 
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where 

H = 2aAo + bBo + cCo - 2dAi . (2.16) 

The Hamiltonian ( p. 161) defines a completely integrable system, since it is effectively one 
dimensional. (Subtract two center elements ( |2.11| ) from four generators and divide by two 
to see that the algebra ( |2.7D - (|2.9| ) is of one degree of freedom.) In the Appendix we show 
how we may eliminate the first term in the Hamiltonian ( |2.16| ) by using the automorphisms 
of the quadratic algebra A. The result is 

H = bBo + cCo - 2dAi , (2.17) 

where the new coefficients b and c are related to o, b, and c through 

a'^ + bc = bc, 2aa + f3b + jc = pi + jc. (2.18) 

Any Hamiltonian that is a function of Ai, Aq, Bq, and Cq will be completely integrable. 
However, for this special case — an arbitrary linear combination — we have formulated the 
problem using the QISM II algebra. It is then possible to apply the method of separation 
of variables, which will be shown in the Section ^ We now give a physical realization of 
the quadratic algebra (^?^)-(|27 



3 The o(4) generalized Lagrange top 



The algebra (2.7)-(p^) may be embedded into the Z//(o(4)) Lie algebra with the six gen- 



erators Jj and Xj, i = 1, 2, 3, and the Poisson brackets 

{e7j,</j} — (^ijkJk ^ {Jii-^j} — ^ijk^k ^ {^ii^j} — ^ijkJk ■ (^-l) 

The lA{o{4)) algebra has the two Casimir elements 

Ci=x'^ + P, C2 = x-J. (3.2) 
The homomorphism between the ^(o(4)) algebra with relations ( |3.1| ) and the rank 1 QISM 



II algebra (2/7)-(^^) is given by the following formulas 



^0 = X1J2 - X2J1 , Bq = -x\ -x\ - JI, Cq = J"^ , Ai = X2,, (3.3) 



where the scalars a, /3, and 7 appearing in (p.4|)-( |2.5| ) are without loss of generahty chosen 
to be 

a = 0, /3 = -l, 7 = 1, (3.4) 

with 6 being equal to 

6 = C2J3 . (3.5) 
Using ( p.lOD -( ^.ll ) we may write the center of the algebra in terms of the o(4) variables 

Qo = Cl + Ci J| , Q2 = Ci + J| . (3.6) 

Notice that from this it follows that J3 commutes with all of Ai, Aq^ Bq, and Cq. It is 
a simple matter to verify equations ( p.7[ )-( [2^ ) from the definitions ( p. 3D and (3.1). The 
integrable Hamiltonian ( [2.17] ) becomes 

H = -bQ2 + {b + c)P + bxl-2dxs, (3.7) 

with J3 as the additional conserved quantity. For the particular values of a, /3, and 
7 in ( |3.4D it can be shown that b and c are always real when a, b, and c are. The 
Hamiltonian ( p. 7] ) corresponds to an o(4) generalized Lagrange top Q. 

In the special case a^+/57 = 0, the algebra A is degenerate. We then have the following 
homomorphism to Z//(e(3)) 

Ao = X1J2 - X2J1 , Bo = -xl-xl, Co = J^, Ai = X3, (3.8) 

where we have put without loss of generality that a = /? = and 7 = 1 with 5 being as 
in ( |3.5D . The algebra generators Jj and Xi satisfy here the Poisson brackets like in ( p.l| ) 
where the last bracket is equal zero (according to the contraction procedure from o(4) to 
e(3)). The first Casimir element Ci = x"^. The Hamiltonian ( |3.7| ) for this case corresponds 
to the Lagrange top in a nonlinear gravity field m. 



4 Separation of variables 



We now return to the Hamiltonian ( ^.16 ) with the Poisson brackets (^.7])- (|2.9| ) and show 



how this system may be integrated by the general method of separation of variables ap- 



plicable to the QISM II algebra. This method was first applied by Sklyanin in |11, |l^] to 
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the representations of the QISM I algebra. See also [g] where the separation of variables 
was applied for the first time to the QISM II algebra under the circumstances close to the 
ones we have in the present paper. 

To do this it is necessary first to apply a similarity transformation to the matrices 
K{u) and T(u) in order to make K{u) triangular. We define the matrix V as 

/ 

V= I , (4.1) 




where a and b are the parameters in K{u) defined by (2.12) and w'^ = + be. Then we 
introduce K{u) and T{u) as 



K{u) = V^K{u){V-y 



T{u) = V-'T{u)V 



A{u) B{u) 
C{u) -A{-u) 



(4.2) 



The matrices K{u) and T[u) are of the same form as K{u) and T{u), respectively. The 
K{u) now looks like 



K(u) 



d 

w ^ — 
u 

c 







d 

w -\ — 
a + w u 



I a + - 

u 

V ^ 



-d + 



d 



(4.3) 



while we have the following relations between the entries of T{u) and T{u) 



a = a + 



a + w 



and 



bBn 



a + w 



(a + w)B[ 







-26a 



Cn 



a + w 



+ (a + w)'j , 



2bAo - b'^Bp + Co 
a + w 



(4.4) 



(4.5) 



where Ai and 6 are not changed. Since det T{u) = det T(n) we have Q2 = Q2 and 
Qo = Qo- Furthermore, K(u) and T{u) satisfy the QISM II algebra too, while t{u) 
defined in ( 2.13| ) is unaltered. Hence, these matrices generate the same integrable system 
as before. Because of the triangular form of K{u) we get 



tiu) = a(i(u) + A{-u)) - -(i(u) - A{-u)) + cC{u) . 



(4.6) 



Now let ui be a zero of C{u), i.e. 



Ciui) = 0., 



(4.7) 



and define 

A+ = -i(-ni), Ar=i(ni). (4.8) 

We may easily evaluate 

X^X^ = -A{ui)A{-ui) = detf{ui) = A{ui) , (4.9) 

where A(u) = detr(n) is given by ( p.lOj ), and from ( |4.6D we get 

t{ui) = w{Xi - Xf) - —{Xi + Xt) . (4.10) 

This is nothing but separation equation for the separation variables ui,Xf. Because f{u) 
satisfies the QISM II algebra we have the important Poisson brackets between the new 
variables ui and (see ^]) 

{m,Xf} = ±Xf, {A+ Ar} = -A'(ni). (4.11) 

From ( [4.g| ) and ( [4.11 ) it follows that we may put 

Xf = ^A{ui) e±Pi , {ui,pi} = l. (4.12) 



Inserting ( [4.12|) into (4.10) and using ( p.l5| ) gives 



H = -2Ja{ui) (w sinhpi + — coshjii^ - {2aa + pb + 7c) uj + 2d-^ . (4.13) 

We have thus formulated the problem as an one-dimensional system with the Hamilto- 
nian ( [4.13| ) in terms of the canonical coordinates ui and pi . The equation of motion looks 
like 

ui = {ui,H}, (4.14) 

from which we have 

{iiif = t{uif - 4A(^.i) (^^ - w^^ . (4.15) 

After the substitution vi = uf the problem can be integrated in terms of elliptic functions. 
The transformation from {Ai, Aq, Bq,Cq) to {ui,pi) is given by equations ( 2.1(]| ), 



and (4.13). To get back we first write A(u) and C{u) in terms of ui and pi through the 



Lagrange interpolation using the data 

C{±ui)=0, A{±ui) = ±Xf, (4.16) 



and the leading parameters a, /3, 7, and 6. This gives 



C{u) = jiu^-uD, (4.17) 



6 uf — u'^ u — ui ^ u + ui 

o ~1~ '^l ~1~ 

u uf 2ui 2ui 



A{u) = - ) + - "\ / + :i__:iA+ + ^^^^Ar ■ (4.18) 



Equating powers of u and using ( 4.12| ) then gives 



5 , v/AR) 



Co = — 7U1 , Aq = —aui — \ A(ui) sinhpi , Ai = ^ H coshpi . (4.19) 

* uf Ul 

To get Bq we insert B{u) = f3v? + Bq into —A{u)A{—u) — B{u)C{u) = A{u) and equate 
powers of u. This gives 



Bq = — \ 2— H 2" + a Ul -\ 2 — sinh pi + 2J A(ni) asinhpi j coshpi 

(4.20) 

Alternatively So may be restored from the Poisson bracket 

{io,ii} = /5C'o-7^o. (4.21) 
We may finally get Ai, Aq, Bq, and Co in terms of ui, pi, Qq, and Q2 by using ( |4.5| ). 



5 Two interacting o(4) tops 

We now generalize the results of the Sections |2| and ^ to the case of not just one but two 
0(4) Lagrange tops interacting with each other. This is done by replacing in ( |2.13| ) K{u) 
with T{u), where T{u) is of the same form as T{u), but with a diff'erent choice of d, /?, 7, 
5, Ai, Aq, Bq, and Cq. Because the QISM II algebra ( |2.lD is closed under the action of 
any similarity transformation, we may write t{u) as 

t{u) = tr f\-u)VT{u)V'^ , (5.1) 



where V is any matrix with determinant 1, while T{u) is given by (p.3|)-( p^ ) and (|3T 
(|3.5D and similarly for T(n). Putting 



V 



a b 
c d 



(5.2) 
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with ad — be = 1 leads to the following integrable system 

t{u) = (a^ + b^ + c^ + -Hu^-G-^, (5.3) 



where 



H = -2{ab + cd)Ao + {a^ + c^)Bo- {b"^ + (f)G^ 





- 2{ac + bd)Ao + {a^ + b^)BQ - (c^ + d^)CQ + 2AiAi , (5.4) 
G = -2(1 + 2bc)AoAo - a'^BoBo - d'^CoCo + 2a{cBoAQ + bAoBo) 

+ b'^CoBo + c^BqCq - 2dibCoAo + cAqCq) + 2{6Ai + 6Ai) . (5.5) 

If we regard the H (^^) as the Hamiltonian then we see that the interaction between the 
two systems is through the term ^ivli. We may use the automorphisms of the quadratic 



algebra A to eliminate the terms in H with Aq and Aq. Using equations ( 2.1^ ) we get the 
following Hamiltonian 

H = A(Co - ^o) + X^HCo - Bo) + 2AiAi , (5.6) 

where A is the single nontrivial parameter in the Hamiltonian satisfying — (A + A^^) = 
+ 6^ + + . The second conserved quantity is 

G = 2{AoAo + 5Ai + 5Ai) - XCqBo - X-^BqCq . (5.7) 

Applying the homomorphism (3.3) gives the final result 



H = j"^ + x + (2;3Cosh6' + X3sinh6')^. (5.8) 

Here A = tanh^, while constant additive terms and multiplicative factors have been ne- 
glected. This is a new integrable system which we call two interacting o(4) Lagrange tops . 
We remark that the number of degrees of freedom is 4 and there are two simple integrals, 



namely: J3 and J3. The fourth integral G (5.7) is of degree 4 in Jj's and Xj's. No separa- 



tion of variables is known for this system or in general for any integrable system given by 
( |5.lD in the case when both T(n) and T{u) are nonscalar. 

The Hamiltonian (|5.8| ) contains interesting subcases. We may easily change the real 
form speaking about two interacting systems on the algebras, for instance, o(2,2)©o(2,2), 
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o(3,l)©o(4) and so on. The further possibihty is the contraction giving the cases hke: 
e(3)©o(4), e(3)©o(3,l), e(3)©e(3) etc. Finally, the two o(4) algebras have two o(3) subal- 
gebras by putting {xi,J2,X3) = (si, 52,53) and similarly (xi, J2,aJ3) = {t2, —tijts). Then 

H = X{sl + 4) + X-\sl + tj) + 2s3t3, (5.9) 
G = X-\siti - Xs2t2f . (5.10) 

This is a special (but non-trivial, one-parameter) case of the general (two parameter) o(4) 
Manakov top g, g. 



6 Further realizations 

In this Section we consider other homomorphisms of the specialized quadratic algebra A 
determined by the condition a^+/37 = 0. Without loss of generality we may put a = (3 = 
and 7 = 1. Applying the techniques of the Section ^ leads to the following system 

H = Co + Co + 2XA1A1 , (6.1) 
G = CoCo + 2X{6Ai + 6Ai-AoAo) + X'^BoBo, (6.2) 

with the single parameter A. We have the following homomorphism of the specialized 
algebra A into Z^(sl(2) © sl(2)) 

^0 = 4{s3t--t3S.), Bo = -m.S-, (6.3) 
Co = (s+ + t+)(s_+t_)-(s3 + t3)2, Ai = 2(t_ -.s_), (6.4) 
5 = 2{s.+t^){Cs-Ct). (6.5) 

Here Cs and Ct are the Casimir elements of the two sl(2) algebras 

Cs = si- s+s- , Ct = tl- t+t- , (6.6) 

and the Poisson brackets for the generators are given by 

{S_,S3} = S_, {S_,S+} = 2S3, {S3,S+} = S+, (6.7) 

{t_,t3} = t-, {t_,t+} = 2t3, {t3,t+} = t+. (6.8) 



11 



We may realize the pair of sl(2) algebras in terms of the canonical variables {x,y,px,Py) 
via the following homomorphism {{x,Px} = 1 etc.) 

xpx x"^ pI 21 

S3 = -—, s. = -, ^^ = Y--2' (6.9) 

where / and m are the values of Cs and Ct, respectively. Inserting these relations into ( |6.l| ) 
gives the following Hamiltonian 

H = \{xpy-ypxf + ^{xpy-ypxf + 2\{x^ -y'^){x^ -y^) 

In polar coordinates x = rcos6, y = rsin^, x = fcos9, y = fs'mO this system becomes 

1 1 ~2 
H = -e^ + -e +2Ar2pcos26'cos26' 



4 4 

I m I fh 



cos^e sin^e cos'^e sin^e' 
G = ^(^ee-2Xr^f'^cos2{e-6) + 2Xr'^f'^cos2{e + e)^ 



(6.12) 



+ 



1-9 / m \ / 1 ~2 I ifi 



4 cos2 6l sin^eiy \4 cos^ sin^ i 
+ 2Xr^f^ ({m - I) cos 26 + {m - I) cos 26 ] 




Regarding H ( |6.12| ) as the Hamiltonian, we can see that we have a system of two particles 
situated in some singular field and interacting to each other through the product of cosines. 
This can be interpreted as two interacting pendulas. This is a new integrable system. The 

and f ^ are constants, being the level values of two simple integrals: s_ + t_ and s_ + 
The fourth integral G ( |6.13| ) is of degree 4 in terms of 6 and 6. We would like to call this 
system the D2 Toda lattice with singular terms. Without the singular terms the system is 
easily integrated by introducing d± = 9±9. However, no integration of the general system 
is yet known. In what follows we will give two more physical examples connected to 
some realizations of the specialized algebra A: namely, the so-called Ai and BC2 rational 



Calogero-Moser systems |14]. 
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The Ai rational Calogero-Moser system is given by the Hamiltonian 

H = lipl + pD + J {p„xj} = 6^j, (6.14) 

with the additional conserved quantity P = pi + P2- Introduce three more variables 

R = xi+X2, Ho = ^{pixi + P2X2) , H+ = ^{xl + xl). (6.15) 

We now have the following realization of the specialized algebra A 

Ao = A{PHo-RH), Bo = 4{P^-4H), Co = A{HH+-HI), (6.16) 
Ax = 2P, (5 = 0, (6.17) 

again with a = /? = and 7 = 1. Since Q2 = 16H, H commutes with all of Aq, Bq, 
Co, and Ai. This observation corresponds to the fact that the Hamiltonian ( |6.14 ) is 
superintegrable p^ . 



Inserting this realization into the Hamiltonian ( p.lj) for two "interacting Ai Calogero- 
Moser systems" leads to the following integrable system 



„2 I ^2 



H = {xiP2 - X2Plf + 2—^— + ixiP2 - P2Xl^^ 

\Xl X2) 

2 , Z.2 



'{xi - X2Y 



+ 2— ^^% + A(pi+p2)(pi+P2). (6.18) 



By the construction of this system, we know that its motion on the levels of two Hamiltoni- 
ans for two A\ Calogero-Moser systems is equivalent to that of the system ( |6.12| ) (because 
both systems give two different realizations of the same specialized algebra A). 

Consider finally the following one-parameter case of the BC2 rational Calogero-Moser 
system given by the Hamiltonian 

1,22^ 1 1 k2 / 1 1 



where k is a constant. This system may be embedded into the specialized algebra A in 
the following manner 

Co = \{Hii^-Bl), Bo = A{N -AH^){N -2H^), (6.20) 
^0 = {iV,Co}, Ai = -2{N -2,H^), 6 = {k^-1)H\ (6.21) 
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where Hq and are given by (|6.15 ) and N is the so-cahed second conserved quantity 
H i{H,N} = 0) 



N 



2 / Pi^pI 



+ 16^^5 — + 



(xf — X^)^ v-^l ~ -^27 •^l-^2 

Because the Hamiltonian H is sitting in 6 it commutes with ah the generators and hence 
the system ( 6.19| ) is sup erintegr able too. In this case the Hamiltonian for "two 

interacting BC2 Calogero-Moser systems" has the following form: 



1 2/ •^\^'2 



(6.22) 



H 



{xiP2 - X2Plf + {xiP2 - X2Plf + 2{xl + xl)U{xi,X2, k) 

+ 2{xl + x^)U{xi,X2, R) + X{N - 3H^){N - 3H^) , 



where 



U{xi,X2,k) 



1 



+ 



(Xl - X2)^ (Xl + X2)^ 



+ 



1 



1 



2 ~^ 2 

1 "^O 



(6.23) 



(6.24) 



It is quite interesting to remark that the motion given by the Hamiltonian ( |6.23| ) on the 
fixed levels of two integrals of motion ( 6.1S| ) (one with tilde and one without tilde) is the 
same as the motion given by the Hamiltonian of the D2 Toda lattice ( |6.12 ). The last two 
results regarding the Calogero-Moser systems were obtained by one of us in p^]. 



7 Two e(3) tops interacting with the Toda lattice 

The next generalization we present is that of inserting the Toda lattice between the two 
e(3) Lagrange tops (see end of the Section 3). This is done by using the monodromy 
matrix for the Toda lattice 

L{u) = L2iu) ■ ■ ■ LN-iiu) , (7.1) 

where 



Li{u) 



^ exp(Q'j 



{qi,Pj} = 5ij. (7.2) 



V 



exp(-gi) u-pi 

The matrix L{u) ( [7.1| ) satisfies the Poisson limit of the QISM I algebra |l^ 

{LW(u),L(2)(t;)} = [r{u - v),L^^\u)L^^\v)] , (7.3) 
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with the same r-matrix (^j) as in (|2.l| ). Now let Ti{u) and T/v(u) be the representations 
described in ( |3.8[) . It is then true that the matrix 

Tiu) = L{u)Tn{u)L-\-u) , (7.4) 



satisfies the QISM II algebra (2J) too |||]. Therefore the trace 



t{u) = tr Tl{-u)L{u)TNiu)L-^{-u) , (7.5) 

is a generating function for the integrals of motion for the An-2 Toda lattice interacting 
with an e(3) Lagrange top at each end. The Hamiltonian for this system is as follows: 

1 ^^'^ 2 2 

^ = o E ^'i + '^i - E exp(gi+i - Qi) - e'^^xi.g + e'^^-^XN,3 ■ (7.6) 

i=2 i=2 

The first three terms describe the kinetic energy of the system, which is a chain of particles 
plus two tops at the ends of the chain. The last three terms are the potentials, the first 
one being of the form of Toda-like interaction between the neighbours in the chain while 
the two last terms describe an interaction of the tops with the chain. It is convenient to 
choose the following representation of the specialized quadratic algebra A. Let a = /3 = 
and 7 = -1. Then (cf. also §, g) 

Ai = coshg, j4o=psinhg, i?o = — sinh^g, (7.7) 



where ci and C2 are arbitrary constants. The Hamiltonian (|7.q ) then becomes 

N N-2 

2 



H = \^Pi - ^ exp(gi+i - Qi) - e^2 cosh qi + e^^-i cosh (7.9) 

i=l i=2 



+ + + + (7 10) 

sinh^ f cosh^ f sinh^ ^ cosh^ ^ ' ^ ' ^ 

This Hamiltonian was proved to be integrable for the first time in p^ ]. It corresponds to 
the most general Toda lattice of the type -D„ with four additional singular terms. Now 
we have given another proof for its integrability and also shown how it appears naturally 
from the combination of the standard An-2 Toda lattice with two tops each interacting 
to the corresponding edge particle of the An-2 Toda lattice. 
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8 Two 0(4) tops interacting with the Heisenberg magnet 

Let us consider the XXX Heisenberg magnet. It is given by the following construction. 
First, we introduce a chain of simple L-operators (A; = 1, . . . , N) 

L,{u) = l + U'' '\] (8.1) 



each of which satisfies the QISM I r-matrix algebra {T{u) = L]t{u)) with the r-matrix 

{r«(n),r(2)(t;)} = [r{u - 7;),rW(n)r(2)(i;)], (8.2) 

where the s-variables have the following Poisson brackets of the sl(2) Lie algebra for any 
k: 

{s^ s^} = ±is^, {s^, s"} = 2is^ (8.3) 

We will also use the real variables s\ and s| defined as = is^. . The L-operators 
( |8.1D have the following properties when conjugate and change the sign of the spectral 
parameter u: 

Lk{-u) = L\.{u) , Lfc(ii) = (j2Lfc(u)fT2 , ueR. (8.4) 

The L-operators (8.1) can be also represented in the form of the scalar product between 
two vectors: 

Lfc(u) = 1 + -(sfc,(T), 4 = (4>Sfc>4)> (8-5) 
u 

where a = (cJi, cr2, CT3) are the standard Pauli matrices. 

The determinants of the L-operators are expressed through the Casimir elements of 
the sl(2) algebras: 

det Lfc(n) = 1 + cl = {slf + {slf + {slf (8.6) 

and have ibzcfc as degeneration points. The monodromy matrix T{u) = Ln{u) ■ ■ ■ Li(n) 
satisfies the algebra (|8.2D too, while its trace 

t{u) =tr T{u), {t{u),t{v)} = 

provides us with the complete set of the integrals of motion in involution for the XXX 
Heisenberg magnet. There is a well-known rule Q to write down the local Hamiltonian for 
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the model. Suppose that we have the homogeneous chain, i.e. all Ck are equal, Ck = c. Then 
one should calculate the logarithm of the generating function t{u) in the common point of 
degeneration of all the L-operators, i.e. u = ic in our case. So, the local Hamiltonian for 
the XXX Heisenberg chain looks like: 

^ 2 

Hioc = log \t{ic)\^ = ^ log [2 + -(4, Sk+i)] . (8.7) 
fc=i ^ 

Let us now proceed further to introducing the boundary conditions for the chain with 
the Hamiltonian ( |8.7[ ). First, we pick up a special representation of the QISM II algebra 
which was introduced in Sections 2 and 3. Suppose we have the Lie algebra so(3)©so(2,l) 
given by the following Poisson brackets for its six generators Sk and tk 

{si,Sj} = EijkSk, (8.8) 

{ti,t2} = -h, {t2,h} = h, {t3,h} = t2. (8.9) 
The Casimir elements are 

Cs = {Slf + {S2f + {ssf = s\ Ct = -ihf - {t2f + {hf = . (8.10) 
The following matrix T{u) 

r(n) = " (8.11) 

V i(n2 + s'^ + f + 253*3 + 2i{siti + 52*2)) 

i{u^ + s2 + ^2 + 253*3 - 2i{siti + 52*2)) \ 
^(53 - t3) - 2.(51*2 - S2h) + ) 

gives a representation of the QISM II algebra ( |2.1| ). 

Let us take the same representation with the tilde-variables and tk which are in 
direct sum to the spins 5 and t: 

{si,Sj} = EijkSk, (8.12) 
{*i,t2} = -t3, {i2,t3} = ii, {t3,ti} = t2. (8.13) 
The Casimir elements are 

Cs = {hf + {§2? + ihf = , Ct = -{hf - {12? + ih? = P . (8.14) 
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The following matrix T(u) 

T{u) = (8.15) 
V i(n2 + s2 + P + 2s3t3 + 2i{siii + hk)) 



i{u'^ + s2 + P + 2s3t3 - 2i{siii + 52^2)) 
u{s3 - h) - 2i(sit2 - hh) + 



gives a representation of the QISM II algebra ( 2.1 



The generating function for the integrals of motion for the XXX Heisenberg chain with 
boundary terms looks like 

t{u) = tr f{u)LN(,u) ■ ■ ■ Li{u)T{u)Li{u) ■ ■ ■ Ln{u). (8.16) 

The local Hamiltonian for the system has the following form: 

N-l 2 

Hio, = logt(ic) = ^ log[2 + -(4,4+i)] (8.17) 
k=i ^ 

+log [C2(S3 - H) - (S3 + h){s^ - + s\{-C^ + S^ + t^ + 2s3t3) 
+2sl{siti + 52*2) + 2sf (Sit2 - S2tl)] 

+log [C^{h - h) - (S3 + h){s^ - P) + s]v(-C^ + S^ + ? + 2s3t3) 
+2s^(siti + 52*2) + 2s^(sif2 - S2tl)] + log . 

This integrable Hamiltonian describes the XXX Heisenberg chain interacting with two 
spins s and t at the one end of the chain and with two spins s and t at the other one. Note 
that in order to get the Hermitean Hamiltonians we had to choose the non-compact so(2,l) 



form of the spins t and i. The boundary terms in the Hamiltonian ( 8.17 ) generalize the 



ones in |15], which described the influence of the external magnetic field and, in quantum 



case, looked like (formula (53) in |jT 



N-l 



H= J2(«+i+'^>n+i+«+i)+b-(^l-b+(^N+C-(yi -c+a^+d^at-d+a+, (8.18) 



n=l 



where b±, c±, and d± are constants. The Hamiltonian like ( ^.18 ) — with six boundary 



terms — was proved recently to be integrable also for the generalization up to the XYZ 
case [ITl. 
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Discussion 



In this paper we constructed several new integrable systems which appear after imposing 
boundary conditions on the known integrable lattices. There were obtained the following 
results: (a) two interacting o(4) tops; (b) a new interpretation of the most general Toda 
lattice of the type (namely, "Dat = An + 2 tops"); (c) quadratic algebra ^ as a 
dynamical algebra of hidden symmetries for the Ai and BC2 Calogero-Moser problems; 
(d) the explicit form of the "local" Hamiltonian ( 8.17] ) for the system which describes an 
interaction of the XXX Heisenberg chain with two o(4,C) tops. 

We believe that our results can be generalized in some possible ways: first, the quan- 
tization which seems quite straightforward, secondly, the g-deformation. There is also an 
open problem to integrate new systems or to separate variables for them. 

It is also very interesting to understand how many spins can be added to the Heisenberg 
chain as some sort of boundary conditons still having the integrability property. In the 
present paper we showed how to add two more spins at each end. Our conjecture is: it 
is possible to extend this up to three spins at each end of the standard spin chain. This 
conjecture came from the study of the hyper geometric orthogonal polynomials |^]. 
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A Automorphisms of the quadratic algebra A 

In this Appendix we discuss the automorphisms of the quadratic algebra A given by the 
relations (|2.7| )-( ^!9[) . It is a property of the QISM II algebra that any similarity 

transformation of a representation T[u) is an automorphism of the QISM II algebra, i.e. 

f{u) =V~^T{u)V , (1.1) 
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for any non-degenerate matrix V satisfies the QISM II algebra too. Rewriting the Ansatz 
(PD-Q as 

r 

T(u) = u^n + Aiul + X + -I, (1.2) 

u 

where I is the 2x2 identity matrix, and 

\ 7 -a ; \Co -Ao J 

it becomes clear that Ai and 5 do not change under the transformation ( |1.1D . We wish to 
find a matrix V that leaves a, P, and 7 unaltered. This will be so if V commutes with 0. 
We consider therefore the following choice 

v{e) = exp i^^n^ , (1.4) 

where = detJl = — (a^ + /37). The only part of T{u) that is changed is X, according 
to the following formula 

x{9) = v{-e)xv{e) . (1.5) 

The relation ( |1.5| ) defines an automorphism of the quadratic algebra A, which leaves a, 
P, 7, 6, and Ai fixed. An arbitrary linear combination of Aq, Bq, and Co may be achieved 
by the following formula 

F = tr{AX) = 2aAo + bBo + cCq , (1.6) 

where 



A 



( 

a c 
[b -a) 



(1.7) 



Substituting X{9) for X in (|L6|) leads to 

F = ti{A{e)X) , (1.8) 

where 

A{9) = V{9)AV{-e) . (1.9) 
This equation thus gives the result of applying the automorphism ( |1.5D to the linear 



combination ( |1.6| ). It is not necessary to know the explicit form of this equation, but 
rather the two main properties 

det A{e) = det A , tr{A{e)n) = tr{An) . (1.10) 
20 



The second relation follows from the fact that V{9) commutes with fi. Letting A(9) = A 
we may write the above relations explicitly as 

+ bc = + bc, 2aa + /3b + jc = 2aa + /Sb + jc. (1.11) 

These relations give two equations to determine the new parameters a, b, and c in terms 
of the old parameters a, b, and c, in agreement with the additional arbitrary parameter 
6. We may thus choose to put any one of the new parameters to zero and thus determine 
the values of the remaining two. 
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